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Braunstein et al. replies: In a comment of the recent paper |T|, Lopez et al. @ 
obtained analytically the short time behavior of the directed percolation depinning (DPD) 
model ||. Their results explains the behavior of the temporal derivative of the interface 
width (DSIW) for all q until a time t ~ e~ 2 . We argue that the fail in reproducing the 
early time regime until the correlations are generated (i ~ 1 at the depinning transition) is 
because the density of active sites of the interface is not a constant p. This density depends 
on time as we will show below. At time t a site i, of a one dimensional lattice of size L, is 
chosen at random with probability 1/L. Let us denote by ht(t) the height of the i-th generic 
site at time t. The set of {hi,i = 1,...,L} defines the interface between "wet" and "dry" 
cells. We shall denote Fi = Fi(hi + 1) the activity of the i-th generic site above the interface. 
If the cell (i,hi + 1) is active Fi = 1 (unblocked) otherwise Fi = (blocked). The time 
evolution for the probability of active sites f(Fi = l,t) = f(t) just above of the interface in 
a time step St = 1/L is 

/(* + «) = f /(*) + (1 " /(*) ■ (1) 

This equation takes into account the probability that a cell above the interface remains 
active after a time step St. The first term take into account the probability of growth f/L 
in the i-th column times the probability p that the new cell of the interface be active. The 
second term is the probability that no growth occurs in the z-th column when the cell is 
active. Taking the limit St — > in Eq. ([I]) and solving the equation, with initial condition 
/(0) = p, we obtain 

f(t)=pe-* t . (2) 

Notice that f(t) is the interface activity density (IAD), i.e. f(t) = {(Fi)} where the brackets 
(braces) denotes averages over the lattice (realizations). Notice that f(t) is close to p until 
t ~ e -2 where the result of Lopez et. al. holds. To obtain a more realistic description for 
the DSIW until the correlation are generated (t ~ 1) it is necessary to take into account 
the temporal dependence of the IAD. Let us consider a growth model in a system of size 
L with density f(t) of active cells in the early time regime. In this regimen the lateral 
correlations are negligible. Assuming independence between Fi and hi, the time evolution 
for the probability of having a column with height h at time t is given by 

p{h,t+st) =p(fc-i,t)M + p(fc,t)(i-M) . (3) 

Lj Lj 

Taking the limit St — > we obtain the master equation for the probability. Using the 
generating function of moments, one can calculate the DSIW: 

£ = /<«). W 

For t •> 1 horizontal correlations are generated and Eq. (|2|) breaks down. From numerical 
simulations we could check that the hypothesis assumed to derive Eq (|3p holds in the neigh- 
borhood of the criticality and in the pinned phase but it breaks down for q q c . However, 
these values of q are not interesting from an experimental point of view HfJ] . A comparison 
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of Eq.(H) with numerical simulations of the DPD model is presented in Fig. [1] for the 
critical value. We can see that until t < e -2 the analytic result of Eq. ([|) and the one 
obtained by Lopez et al. are coincident with the numerical results of the DSIW. This is 
because in this regime f(t) ~ p. As time goes on, f(t) decays and the hypothesis of Lopez 
et al. does not hold as we can see in this Figure. However, our analytical result predicts the 
DSIW until t~l. 

Finally, we argue that contrary to what it is claimed in the comment of Braunstein 
and Buceta's formula does describe the macroscopic behavior of the interface even when the 
solution is not exactly a matching between the early [Eq. (|^)] and the asymptotic regimes. 
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FIGURES 

FIG. 1. Numerical results for the DSIW (circles) and the IAD (triangles) for the DPD model 
in a system of size L = 8192 for q c = 0.539. Continuous lines correspond to Eq. (|j) and fit the 
data until t ~ 1. Dashed line correspond to the approximation of Lopez et ah. 
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